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ABSTRACT 


A  caustic  flooding  study,  using  Ottawa  sand  and  native 
crude  oil,  was  conducted  to  investigate  the  effects  of 
linear  rate  and  adsorption  on  oil  recovery.  Thirteen  runs 
were  conducted,  twelve  with  the  sodium  hydroxide  solutions 
displacing  the  Wainwright  crude  oil  and  a  final  displacement 
of  an  artificial  brine  by  a  0 . 1  %  by  weight  sodium  hydroxide 
solution.  All  displacements  were  carried  out  at  23°C. 

Three  different  sodium  hydroxide  concentrations  were 
used  and  an  improvement  in  oil  recovery,  over  that  obtained 
for  a  brine  flood,  was  noted  for  the  0.01  and  0.1  %  by 
weight  displacements  but  not  for  the  2.0  %  displacement.  Oil 
recovery  at  two  pore  volumes  injected  was  superior  in  the 
0.01  %  displacement,  while  the  breakthrough  recovery  was 
superior  in  the  0.1  %  displacements.  Linear  rate  had 
differing  effects  on  the  0.01  and  0.1  %  displacements. 

The  oil  recovery  in  the  0.1  %  by  weight  displacements, 
where  the  recovery  mechanism  appeared  to  be  emulsification 
and  entrainment,  seemed  essentially  independent  of  rate. 

Only  at  the  extreme  rates,  0.03  and  14.49  m/day,  was  the 
recovery  significantly  affected. 

The  oil  recovery  in  the  0.01  %  by  weight  displacements, 
where  the  recovery  mechanism  was  postulated  to  be 
emulsification  and  entrapment,  exhibited  a  maximum  of  75.5  % 
of  the  initial  oil  in  place  at  3.05  m/day. 

Consumption  of  sodium  hydroxide  by  the  sandpack  was 


found  to  be  significant. 
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1.  INTRODUCTION 


Large  quantities  of  viscous  crude  oil  exist  throughout  the 
world;  maximizing  the  recovery  from  these  deposits  is  of 
great  interest  in  these  days  of  ever - i ncreas i ng  oil  prices. 
It  is  well  Known  that  conventional  water  flooding,  although 
increasing  the  recovery  of  oil,  leaves  a  substantial 
residual  oil  saturation.  The  possibility  of  reducing  this 
residual  saturation  via  the  use  of  chemical  additives  is,  as 
a  result,  being  actively  investigated. 

Caustic  is  just  one  of  the  many  possible  additives  that 
has  been,  and  continues  to  be,  studied.  Although  many 
variables,  such  as  interfacial  tension,  are  Known  to  affect 
the  recovery  process,  little  theoretical  worK  has  been 
attempted,  and  no  satisfactory  mathematical  representation 
of  the  process  exists  at  this  time. 

The  present  study  was  devoted  to  the  study  of  the 
effect  of  rate  on  oil  recovery  in  caustic  floods,  and  the 
adsorption  of  caustic  on  the  rocK  surfaces.  It  is  hoped  that 
a  better  understanding  of  the  adsorption  phenomenon  could 
lead  to  a  more  complete  theoretical  frameworK  for  caustic 
f 1 ood i ng . 


1.1  Organization  of  the  Thesis 

A  review  of  the  major  worKs  relating  to  caustic 
flooding  is  presented  in  Chapter  2.  The  vast  quantity  of 
literature  on  the  many  aspects  of  caustic  flooding 
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necessitates  that  this  is  not  a  comprehensive  review,  but 
rather  a  review  of  selected  works.  The  proposed  recovery 
mechanisms  are  presented,  and  some  of  the  factors  affecting 
them  examined. 

Some  of  the  factors  effecting  caustic  flooding  are 
given  in  the  third  section  of  this  chapter.  The  recently 
reported  possibility  of  time  dependence  is  briefly  examined 
in  the  fourth  section. 

The  final  section  of  this  chapter  is  devoted  to  the 
basic  principles  of  adsorption  and  dispersion  theory  as 
applied  to  miscible  displacements.  The  evaluation  of  the 
various  rate  and  dispersion  constants  is  examined. 

The  third  chapter  states  the  various  specific  areas 
that  the  study  investigated. 

The  materials  and  experimental  technique  employed 
during  this  study  are  outlined  in  the  fourth  chapter.  As  the 
techniques  used  are  routine,  they  are  mentioned  only 
briefly. 

Presented  in  Chapter  5  are  the  results  of  the  tests. 
These  results  are  discussed  in  an  attempt  to  relate  them  to 
the  currently  accepted  mechanisms  of  caustic  flooding.  The 
results  and  conclusions  of  this  study  are  summarized  in 
Chapter  6. 

The  final  chapter  outlines  specific  recommendations  for 
further  study.  It  is  hoped  that  work  in  these  areas  will 
provide  further  insight  into  the  mechanisms  of  caustic 
f 1 ood i ng . 


2.  THEORY  AND  LITERATURE  REVIEW 


2.1  History  of  Chemical  Water  Flooding 

The  use  of  chemical  additives  in  flood  water  was 
proposed  by  A tki nson ( 1 927 ) .  The  benefits  of  alkaline  flood 
water  additives  were  known  even  before  this  as  shown  by  the 
works  of  Nut t i ng ( 1 925 ) ,  Uren  and  Fahmy(1927)  and  Bekstrom 
and  Van  Tuy 1 ( 1 927 ) . 

Nut t i ng ( 1 925 )  described  the  addition  of  alkali  salts, 
such  as  sodium  carbonate,  to  the  flood  water  in  order  to 
improve  waterflood  recovery.  He  dismissed  the  use  of 
stronger  bases,  sodium  hydroxide  for  example,  on  the  grounds 
that  they  would  be  excessively  reactive  with  the  crude  oil 
and  would  be  used  up  by  the  oil  before  they  could  be 
effective.  Despite  this,  Atkinson's  patent  of  1927  described 
the  benefits  of  stronger  bases  including  sodium  and 
potassium  hydroxide.  Bekstrom  and  Van  T uy 1 ( 1927)  also 
reported  improved  oil  recovery  with  the  use  of  both  strong 
and  weak  bases.  Uren  and  Fahmy(1927)  confirmed  the  results 
of  Bekstrom  and  Van  T uy 1 ( 1927)  and  attributed  the  increase 
to  the  release  of  trapped  oil  from  the  rock  surface. 

Nut t i ng ( 1 925 )  had  found  that  the  alkali  prevented  the 
formation  of  semi-solid,  crude-oi 1 -water  interfacial  films; 
however,  he  did  not  correlate  this  with  increased  recovery. 
Atki nson ( 1 927 ) ,  on  the  other  hand,  recognized  that  the 
capillary  and  adhesive  properties  of  the  oil  would  cause  it 
to  be  held  in  the  pore  structure,  and  that  the  addition  of 
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an  alkali  overcame  these  forces  and  released  the  oil.  This 
mechanism  was  postulated  to  be  a  combination  of  wettability 
alteration  and  interfacial  tension  reduction  in  agreement 
with  the  findings  of  Bekstrom  and  Van  Tuyl(1927). 

Bekstrom  and  Van  T uy 1 ( 1927)  also  maintained  that  the 
formation  of  emulsions  was  undesirable.  This  conflicts  with 
the  views  of  Subkow( 1 942 ) ,  who  considered  emulsification  to 
be  the  essential  first  step  in  the  recovery  process.  He 
stated  that  emulsification  must  be  followed  by  the 
entrainment  of  the  emulsion  in  the  flowing  alkali  stream  in 
order  to  increase  oil  recovery. 


2.2  Proposed  Recovery  Mechanisms 

Johnson ( 1 976 )  performed  a  thorough  review  of  the  status 
of  caustic  flooding  methods  and  classified  the  proposed 
mechanisms  into  the  following  four  types: 

1)  Emulsification  and  Entrainment 

2)  Wettability  Reversal  (oil-wet  to  water-wet) 

3)  Wettability  Reversal  (water-wet  to  oil-wet) 

4)  Emulsification  and  Entrapment 
Each  of  these  mechanisms  is  reviewed  below. 

2.2.1  Emulsification  and  Entrainment 

This  mechanism  was  first  proposed  by  Subkow(1942)  in 
relation  to  the  recovery  of  bitumen.  He  stated  that  the 
formation  of  an  emulsion  was  the  essential  first  step  in  the 
recovery  process,  followed  by  the  entrainment  of  this 
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emulsion  in  the  flowing  alkali  with  the  mixture  subsequently 
produced.  He  postulated  that  the  alkali  solution  caused  the 
formation  of  emulsification  agents  which  lowered  the 
interfacial  tension.  These  agents  were  formed  from  reactions 
between  the  sodium  hydroxide  and  the  organic  acids  present 
in  the  bitumen.  The  bitumen  was  emulsified  into  an 
oil-in-water  emulsion  and  stripped  from  the  sand  grains.  He 
foresaw  the  dangers  of  using  high  concentrations  of  alkali 
on  two  grounds: 

1)  Excessive  alkali  could  cause  the  formation  of  a 
water-in-oil  emulsion. 

2)  The  high  electrolytic  concentrations  could 
prevent  the  formation  of  the  emulsion. 

Reisberg  and  Doscher ( 1 956 )  studied  the  effect  of  pH  on 
the  oil-water  interface  using  a  Ventura,  California,  crude 
oil.  They  concluded  that  for  the  same  interfacial  tension  to 
oil,  a  caustic  solution  would  displace  the  oil,  while  an 
acid  solution  would  not.  They  attributed  this  to  the 
stability  of  the  rigid  films  at  the  interface  confirming  the 
earlier  work  of  Nut t i ng ( 1 925 ) .  The  rigid  film  was  isolated 
and  found  to  be  a  highly  oxygenated,  low  molecular  weight 
constituent  of  the  resins  and  asphaltenes  found  in  the  oil. 
This  confirmed  the  earlier  work  of  Dodd,  et  a  1 . ( 1952) . 

Reisberg  and  Doscher ( 1 956 )  felt  that  the  ability  of  the 
caustic  to  prevent  adherence  of  the  oil  to  the  rock  surfaces 
was  important;  however,  they  placed  primary  importance  on 
the  emulsification  and  entrainment  mechanism.  They  concluded 
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that  caustic  flooding  was  impractical  because  adsorption  on, 
and  reaction  with,  the  rock  and  displacement  of  connate 
water  would  cause  the  alkali  to  fall  behind  the  oil -water 
interface.  As  a  result,  they  concluded  that  an  increase  in 
oil  recovery  would  not  be  evident  until  several  pore  volumes 
of  caustic  had  been  injected,  an  economically  prohibitive 
situation. 

Chant  1979),  working  with  Wainwright  crude  oil, 
concluded  emulsification  of  the  oil  to  be  vital  for 
increased  oil  recovery.  This  conclusion  was  based  on  the 
insignificant  recovery  increase  in  tests  where 
emulsification  did  not  occur. 

2.2.2  Wettability  Reversal  (Oil-wet  to  Water-wet) 

Wagner  and  Leach(1958)  obtained  improved  oil  recovery 
through  the  injection  of  chemical  solutions  which  reversed 
rock  wettability  from  oil -wet  to  water-wet.  The  chemicals 
included  acids,  bases  and  certain  salts.  They  recognized 
that  the  improvement  could  be  attributed  to  the  improved 
mobility  ratio,  due  to  favourable  changes  in  oil  and  water 
relative  permeabilities,  in  the  region  where  oil  was  still 
flowing.  This  decrease  in  mobility  ratio  could  also 
temporarily  halt  the  gradual  increase  in  producing  water-oil 
ratio. 

They  reasoned  that  the  injected  water  would  displace 
the  connate  water  in  front  of  it  and,  thus,  the  injected 
chemical  would  only  contact  the  residual  oil  left  behind  the 
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connate  water  bank.  Since,  in  a  water-wet  system,  the 
residual  oil  is  discontinuous  they  concluded  that  a 
water-wet  system  would  not  respond  to  this  wettability 
change . 

The  acidic  solutions  were  dismissed  on  the  grounds  that 
they  are  too  reactive  with  most  reservoir  rocks. 

Consequently  Leach,  et  a  1  . ( 1962)  performed  both  laboratory 
and  field  tests  using  sodium  hydroxide  and  obtained 
improvements  in  oil  recovery  in  the  regions  where  the 
wettability  change  took  place.  These  improvements  were 
similar  to  those  obtained  by  Wagner  and  Leach(1958). 

Mungan ( 1 966 . a )  performed  laboratory  tests  similar  to 
those  of  Leach,  et  a  1 .  (  1962) ,  and  confirmed  the  improvement 
in  mobility  ratio.  He  also  found  that,  for  his  system,  the 
process  was  temperature  dependent.  The  recovery  mechanism 
worked  well  at  7 1 . 1° C  ( 1 6  0  °  F )  but  not  at  21. TC  (70°F).  This 
temperature  dependence  was  later  confirmed  by  Cooper ( 1 97 1 ) , 
and  attributed  to  the  temperature  sensitivity  of  interfacial 
tension.  A  later  paper  by  Mungan ( 1 966 . b )  confirmed  the 
benefits  of  a  wettability  alteration  from  oil-wet  to 
water-wet  on  the  recovery  of  oil  from  Teflon  cores. 

Ehrlich,  et  a  1 . ( 1974)  concluded  that  wettability 
reversal  was  a  viable  process  in  light  oil  reservoirs.  They 
noted,  however,  that  if  the  interfacial  tension  was  too  low, 
oil  recovery  would  still  increase,  but  the  mechanism  would 
be  emulsification  and  entrainment. 

Scott,  et  a  1 . ( 1965)  found  a  significant  increase  in  oil 
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recovery  with  the  use  of  weak  caustic  solutions.  They 
attributed  the  increase  in  oil  recovery  to  a  shift  in 
wettability  toward  water-wet.  Their  core  was  initially 
"neutral  or  slighty  oil-wet".  The  observed  shift  in  the 
unsteady  state  relative  permeability  curve,  toward  increased 
oil  permeability,  was  offered,  as  evidence  for  the  proposed 
mechani sm. 

Minssieux( 1977 )  presented  similar  evidence  to  support 
his  conclusion  that  the  addition  of  caustic  increased  the 
wettability  to  water  of  the  reservoir  rock  surfaces.  He  did 
not,  however,  report  the  original  state  of  his  system. 

2.2.3  Wettability  Reversal  (Water-wet  to  oil-wet) 

Cooke,  et  a  1 . ( 1974)  reported  a  third  mechanism  by  which 
sodium  hydroxide  could  improve  oil  recovery.  They  showed 
that,  under  certain  conditions  of  pH,  salinity  and 
temperature,  some  systems  could  be  converted  from  water -wet 
to  oil -wet.  With  the  right  crude  oil,  this  reversal  would  be 
accompanied  by  a  dramatic  reduction  in  interfacial  tension. 
The  low  values  of  interfacial  tension  are  due  to  the 
reaction  of  the  alkali  with  organic  acids,  at  the  water-oil 
interface,  forming  soaps.  The  solubility  of  these  soaps,  in 
the  aqueous  phase,  is  highly  dependent  on  its  salinity.  The 
steps  in  the  process  are: 

1)  Conversion  of  the  rock  from  water-wet  to 
oi 1 -wet . 

2)  The  previously  discontinuous  residual  oil 
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saturation  becomes  continuous,  and  oil  is  able 
to  flow. 

3)  The  low  interfacial  tension  induces  the 
formation  of  oil-in-water  emulsions. 

4)  The  emulsion  droplets  block  off  the  smaller  pore 
throats  causing  large  pressure  drops  in  the 
region . 

5)  The  high  pressure  gradients  overcome  the 
capillary  forces,  already  lessened  by  the  lower 
interfacial  tension.  This  further  reduces  the 
residual  oil  saturation. 

The  result  is  that  a  high  water  content  emulsion  is  left 
behind  in  the  pore  structure. 

The  most  distinctive  feature  of  this  mechanism  is  the 
wettability  change.  In  this  regard,  as  can  be  seen  from 
Figure  1,  the  pH  and  salinity  of  the  aqueous  phase  appear 
critical . 

The  work  of  Mi nss i eux ( 1 977 ) ,  however,  showed  that  when 
salt  was  present  in  the  caustic  solution  the  interfacial 
tension  was  not  as  low  as  that  obtained  in  its  absence.  The 
Emeraude  crude  used  in  part  of  his  study  typified  this 
finding;  however,  the  improvements  obtained  in  oil  recovery 
were  still  attributed  to  a  wettability  shift  toward  oil-wet. 
This  would  suggest  that  only  part  of  the  proposed  mechanism 
is  required  to  improve  oil  recovery. 
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FIG  1  :  EFFECT  OF  FLOOD  WATER 
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2.2.4  Emulsification  and  Entrapment 


This  final  mechanism  was  proposed  by 
Jennings,  et  a  1 . ( 1974) ,  and  is  restricted  to  water-wet 
systems.  In  this  mechnism  the  emulsion  is  not  produced  as  in 
emulsification  and  entrainment,  rather,  it  blocks  off  the 
smaller  pore  throats  and  increases  vertical  and  areal  sweep 
efficiencies.  The  residual  oil  saturation  is  not 
substantially  reduced  and,  as  a  result,  this  mechanism  is 
best  suited  to  viscous  oils  where  sweep  efficiency  is 
norma  1 1 y  poor . 

This  mechanisim  involves  the  residual  oil  being 
emulsified  in  situ,  due  to  the  low  interfacial  tension,  and 
being  carried  by  the  flowing  alkali  until  the  oil  droplets 
are  trapped  in  the  pore  throats.  This  results  in  reduced 
water  mobility  that  improves  both  vertical  and  areal  sweep 
efficiencies.  Due  to  the  fact  that  the  emulsified  oil 
droplets  are  quickly  trapped,  the  residual  oil  saturation  is 
not  reduced. 

The  two  possible  reasons  the  oil  droplets  are  not 
pushed  through  the  throats  are  : 

1)  The  interfacial  tension  is  not  low  enough  and/or 

2)  The  pressure  gradient  is  not  high  enough. 

2.3  Factors  Effecting  Caustic  Flooding 

Various  factors  have  been  reported  to  effect  the 
caustic  flooding  process.  These  include: 

1)  Crude  Oil  Composition:  It  is  agreed  that  the 
nature  of  the  polar  compounds  in  the  oil 
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determines  whether  the  caustic  will  aid  the 
recovery  process.  A  good  review  of  the  subject 
has  been  published  by  Ehr 1 i ch ( 1 974 ) . 

Mi nss i eux ( 1 977 )  suggested  that  the  acid  index 
should  be  greater  than  about  0.5  mg  KOH/g  of 
crude . 

2)  Water  Composition:  This  is  also  an  important 
variable  as  significant  amounts  of  multivalent 
positive  ions,  such  as  calcium,  can  cause 
excessive  consumption  of  the  caustic. 

3)  Rock  React i vi ty :  Minssieux( 1977 )  suggested  that 
the  mineralogy  of  the  rock  does  not  impose  any 
special  limitations,  except  in  the  case  of  large 
amounts  of  gypsum.  Johnson ( 1 976 ) ,  on  the  other 
hand,  states  that  the  rock  is  important  as  it 
accounts  for  both  consumption  and  adsorption  of 
the  caustic,  in  addition  to  playing  a  major  role 
in  deciding  which  recovery  mechanism  will 

domi nate . 

4)  Caustic  Concentration:  It  appears  that  the 
optimal  level  of  caustic  concentration  depends 
upon  the  prevailing  recovery  mechanism. 

Johnson ( 1 976 )  reports  concentrations  as  low  as 
0.001  weight  percent  for  emulsification 
mechanisms  and  as  high  as  15  weight  percent  for 
wettability  reversals. 
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2.4  Time  Dependence  of  Caustic  Flooding 

Recent  work  by  Chan(1979)  has  suggested  the  possible 
need  of  a  minimum  residence  time  for  the  caustic  to  react. 

He  suggested  that  the  formation  of  the  oil -in-water 
emulsions  was  dependent  on  the  time  during  which  the  caustic 
was  in  contact  with  the  system.  He  offered  no  explanation 
for  this. 

Minssieux( 1977 )  studied  the  effect  of  time  on  the 
interfacial  tension  between  caustic  solutions  and  various 
crude  oils.  He  found,  with  a  Lacq  crude  (acid  index  0.5  mg 
KOH/g  crude),  that  the  interfacial  tension  increased  with 
time  for  any  particular  caustic  concentration.  An  Emeraude 
crude  (Acid  index  1 . 2  mg  KOH/g  crude)  showed  the  same  trend. 

Reisberg  and  Doscher ( 1 956 )  also  studied  the  effect  of 
time  on  interfacial  tension  against  deionized  water  for  5 
crude  oils  and  purified  benzene.  In  all  cases,  they  reported 
a  rapid  decrease  in  interfacial  tension  followed  by  a 
gradual  reduction  toward  some  equilibrium  value. 

McCaf fery ( 1 975 )  studied  interfacial  tension  and  aging 
behaviour  of  4  crude  oils.  He  found  that  the  interfacial 
tension  between  the  oils  and  sodium  hydroxide  solutions 
increased  with  the  interface  age.  The  oils  he  studied 
included  Viking  Kinsella  Wainwright  B,  which  displayed  an 
interfacial  tension  minimum  at  about  0.25  wt .  %  sodium 
hydroxide  in  injection  water.  This  injection  water  contained 
a  total  calcium  and  magnesium  ion  content  of  approximately 
12  grams/m3  (12  parts  per  million).  For  this  system,  the 


' 


14 


interfacial  tension  increased  from  about  8x10'5  N/m 
(0.08  dynes/cm)  to  about  4x10-3  N/m  (4  dynes/cm)  over  some 
70  minutes.  An  aqueous  phase  of  limestone  equilibrated 
water,  containing  some  26  grams/m3  (26  parts  per  million) 
calcium  and  magnesium  ions,  displayed  a  similar  initial 
interfacial  tension,  but  increased  more  rapidly  with  aging, 
reaching  1 x 1 0  ~ 2  N/m  (10  dynes/cm)  after  only  20  minutes. 

McCaf fery ( 1 975 )  also  observed  that  the  aging  period 
increased  with  the  size  of  the  oil  droplet  used.  This  was 
believed  to  account  for  the  longer  duration  of  a  low 
interfacial  tension  for  the  injection  water  system. 

It  was  found  that  the  sodium  hydroxide  in  the  aqueous 
phase  was  not  noticeably  consumed  during  the  aging  periods. 


2.5  Dispersion  and  Adsorption 

2.5.1  Introduct ion 

Miscible  displacement  has  been  the  subject  of  intensive 
investigation  due  to  its  great  potential  in  tertiary  oil 
recovery.  Although  caustic  flooding  is  not  a  miscible 
process,  the  theory  developed  for  miscible  and  mi sci ble- type 
displacements  allows  the  study  of  the  adsorption  phenomenon. 
Adsorption  is  believed  to  play  a  major  role  in  the  caustic 
recovery  mechanisms  but  little  theoretical  work  on  the 
adsorption  of  caustic  on  rock  surfaces  has  been  performed.  A 
study  of  a  miscible  caustic  flood,  (e.g.  caustic  displacing 
brine)  should  provide  much  needed  insight  into  the 
phenomenon . 
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2.5.2  Dispersion  Theory 

In  any  miscible  displacement  process,  the  mixing  of  two 
fluids  is  influenced  by  diffusion  and  dispersion  phenomena. 
Dispersion  is  caused  by  the  uneven  movement  of  fluid,  or  by 
the  concentration  gradient  between  the  two  fluids.  Two  types 
of  dispersion  may  take  place  in  a  porous  medium;  the  first 
is  longitudinal,  which  is  in  the  direction  of  flow,  and  the 
second  is  transverse,  which  is  perpendicul ar  to  the 
direction  of  flow.  The  "diffusion  equation"  for  longitudinal 
dispersion  is  well  understood  and  has  been  solved  by  several 
workers,  including  Brigham,  et  a  1 . (  1 9 6  1  )  ,  using  appropriate 
boundary  conditions.  One  of  the  solutions  presented  is: 

xi 

1  -  erf  -  (  1  ) 

2/Dt_ 

where : 

x-j  =  distance  from  the  midpoint  of  the  flood  front,  m 
C  =  concentration  of  solute  in  the  fluid  phase,  kg/m3 
Cg  =  concentration  of  solute  in  input,  kg/m3 
D  =  longitudinal  dispersion  coefficient,  m2/sec 
t  =  time,  seconds 

As  reported  by  Brigham,  et  a  1 . ( 1 95 1 ) ,  there  is  a  simple 
relationship  between  x-j  and  the  distance  from  the  inlet  of 
the  porous  medium  x,  as  shown  below: 


x  -  ut 


x  -  (L/T)t 


(2) 
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where : 

L  =  length  of  the  porous  medium,  m 
T  =  time  to  inject  one  pore  volume,  seconds 
u  =  average  pore  velocity,  m/sec 
Equation  (1)  gives  the  well  Known  "S"  shape 
concentration  profile  in  the  porous  medium  at  a  given  time. 
Experimentally,  concentration  data  is  obtained  at  a  single 
point,  normally  the  outlet,  therefore,  some  modification  of 
Equation  (1)  is  needed  to  make  it  convenient  for 
experimental  purposes.  Brigham,  et  a  1 . ( 1961 )  have  suggested 
the  fol lowing : 

V  -  V 

x,  =  L  -  (  3 ) 

p 

t  =  if  (4) 

P 


where : 

V  =  total  pore  volume  of  sand  pack,  ml 

r 

V  =  volume  of  liquid  recovered  at  time  of  sample,  ml 
Then  the  argument  of  the  error  function  may  be  rewritten  as, 


=  L^Vp  - 
2/Dt  2 /DOT 


2/DlT 


U 


(5) 


where : 


U  = 


VP  -- 

/V 
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The  error  function  parameter,  U,  may  be  used  to  predict  the 
growth  of  the  transition  zone.  A  plot  of  U  versus 
concentration  on  arithmetic  probability  paper  should  result 
in  a  straight  line. 

Brigham,  et  a  1 . ( 1 96 1 )  obtained  the  following 
relationship  between  the  dispersion  coefficient,  D,  and  the 
error  function  parameter  U: 


D  = 


1 


V  T 
P 


U90  ~  U10 

3.625 


(6) 


where : 

UgQ  =  Value  of  U  at  C/Cq  =  0.9 
U-jq  =  Value  of  U  at  C/CQ  =  0.1 
Equation  (6)  allows  evaluation  of  D  from  experimental  data. 

Raimondi,  et  a  1 . ( 1959)  presented  the  following 
relationship  between  the  longitudinal  dispersion  coefficient 
D,  and  the  mixing  coefficient  a  : 


D  =  D*  +  aU  ( 7 ) 

The  mixing  coefficient,  a  is  dependent  upon  core  lithology 
and  is  independent  of  fluid  properties.  At  high  flow  rates, 
the  parameter  D*  ,  the  diffusion  coefficient,  may  be 
neglected.  Kasraie( 1979 )  found  that  the  mixing  coefficient 
was  rate  sensitive  with  a  maximum  occurring. 
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2.5.3  Simultaneous  Dispersion  and  Adsorption 

Banks  and  A 1 i ( 1964)  modified  the  dispersion  equation  to 
include  the  adsorption  phenomenon.  The  resulting  mass 
balance  on  an  incremental  volume  yields: 


aC  ^  1  -  <j>  9Cs  .  aC  _  n  82C 
3t  *  at  ax  3x2 


(8) 


The  general  form  of  the  rate  equation  expressing  solute 
transfer  between  the  solid  and  fluid  phases  is: 


aC 

at1  ■  f<c*Cs> 


(9) 


Once  the  precise  form  of  the  rate  equation  is  specified, 
equations  (8)  and  (9)  may  be  solved  for  specified  initial 
and  boundary  conditions. 

Rate  equations  commonly  used  in  this  solution  are: 


(10) 


and 


C  =  IOC  +  K',  Ml) 

S  I  c 

where  K'^,  K^,  K-j  and  K2  are  rate  constants.  Equation  (11) 
indicates  the  condition  of  local  equilibrium  between  the 
fluid  and  solid  phases  at  all  points  in  the  porous  medium. 
This  is  a  special  case  of  Equation  (10)  which  is  the 
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nonequilibrium  relationship.  Carslaw  and  Jaeger(1959)  solved 
the  equilibrium  case  for  the  following  boundary  and  initial 
condi t ions : 


C(x,  0)  = 

0 

(  12) 

C$(x,  0)  = 

0 

(13) 

C(0,  t)  = 

co 

(  14) 

The  solution  of  Equation  (7) 

wi th  Equat ions  (11), 

( 12) ,  (13) 

and  (14)  is: 


=  *(£,  n) 

0 


where : 


$(?,  n) 


1 

2 


erfc 


1  - 


^2/^n 


exp 


and 


5 


n 


e 


ut 

ex 


_D_ 

ux 


K 


1 


(15) 


(16) 


(17) 

(18) 

(19) 
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The  general  solution  of  the  nonequilibrium  case  is 
obtained  from  Equations  (8),  (10),  (12),  (13)  and  (14).  For 


this  case  Lapidus  and  Amundson ( 1 952 )  obtained: 


C 


0 


exp 


dt 


(20) 


0 


in  which 


x 

4DA  " 


G(t) 


e 


e 


dA  (21) 


where 


(22) 


and  Iq  is  the  modified  Bessel  function  of  the  first  Kind. 

Banks  and  A 1 i ( 1964)  presented  two  special  cases  of  this 
general  solution.  The  first  assumes  K-j  *  0 ,  =  0  and 

D  t  0.  With  these  assumptions,  the  solution  reduces  to  the 
form: 


+ 


(23) 


in  which 


ut 

X 


(24) 


. 


21 


n 


(25) 


and 


o 


1  +  l  ~ 

4> 


u 


(26) 


The  second  case  assumes  *  0 ,  *  0  and  D  =  0.  In 

this  instance,  with  adsorption  but  no  dispersion,  Banks  and 
A 1 i ( 1964)  report  the  result  to  be: 


-  J(03,  Ip ) 

Lo 


in  which 


J(o),  l|>) 


and 


(0 


K-j  x 
u 


' P  ~  ^2  u"  (^o  "  ^ 


(27) 


(23) 


(29) 

(30) 


2.5.4  Application  to  Experimental  Results 

The  results  presented  in  the  previous  two  sections 
provide  a  theoretical  framework  within  which  to  examine 
experimental  data.  The  major  difficulties  in  attempting  to 
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match  the  theory  to  the  results  lie  in  the  evaluation  of  the 
rate  constants  and  dispersion  coefficient.  Assigning  values 
to  these  quantities  involves  the  use  of  history  matching 
techniques,  and  consequently  invokes  the  assumptions 
embedded  in  the  theories.  That  is,  it  is  required  to  assume 
the  validity  of  the  theory  prior  to  any  attempt  to  match  it 
with  experimental  data. 

2.5.4.  1  Evaluation  of  the  Dispersion  Coefficient 

The  dispersion  coefficient,  D,  may  be  easily  obtained 
using  the  procedure  of  Brigham,  et  a  1 . ( 1 96 1 ) ,  outlined  in 
Section  2.5.2.  An  experimental  system  where  adsorption  is 
not  present  must  be  used,  as  the  theory  upon  which  their 
method  is  based  is  for  dispersion  only.  Since  the  dispersion 
coefficient  is  dependent  solely  on  core  properties,  a  value 
obtained  from  a  simple  dispersion  system  may  subsequently  be 
used  in  a  dispersion  and  adsorption  system.  The  same  flow 
rate  must  be  used  in  both  experiments  since  Equation  (7) 
shows  that  the  relative  effects  of  diffusion  and  dispersion 
are  velocity  dependent. 

2.5. 4. 2  Evaluation  of  the  Rate  Constants  for  the 

Nonequilibrium  Case 

The  values  of  the  nonequi  Ibr  i  urn  rate  constants,  K-j  and 
,  may  both  be  evaluated  from  the  experimental  data  via 
the  second  special  case  of  Banks  and  A 1  i (  1964) ,  presented  in 
Section  2.5.3.  The  assumption  for  this  case',  that  D  =  0, 
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does  not  affect  the  constant  evaluation,  as  might  first  be 
suspected.  In  the  absence  of  adsorption,  it  is  clear  that, 
with  or  without  dispersion,  the  concentration  distribution 
curve  must  pass  through  the  point  C/CQ  =  0.5,  Co  =  1.0. 
Therefore  any  delay  in  reaching  this  point  is  due  entirely 
to  adsorption  effects. 

K-j  may  be  evaluated  by  point  matching  at  Co  =  1.0;  from 
Equations  (27)  and  (30)  we  obtain: 


J(u,  0) 


(31  ) 


In  this  special  case,  the  J-function  reduces  to  a  simple 
exponential,  and  Equation  (31)  becomes: 


Thus,  K-j  is  simply  evaluated  from  the  above  equation.  Having 
evaluated  K-j  ,  K2  may  be  evaluated  by  point  matching  at 
C/Cq  =  0.5.  This  may  be  done  either  graphically  (e.g.  from 
Figure  2  from  Banks  and  A 1 i ( 1 964 ) )  or  numerically. 
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3.  STATEMENT  OF  THE  PROBLEM 


The 

of  rate 
med i um . 

1  . 

2. 
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objective  of  this  study  was  to  examine  the  effect 
on  caustic  displacements  in  an  unconsolidated  porous 
Specifically,  it  was  desired  to  investigate: 

The  recovery  of  Wainwright  crude  oil  as  a  function 
of  injection  rate  for  a  continuous  sodium  hydroxide 
di spl acement . 

The  sensitivity  of  the  rate  dependence  to  the 
concentration  of  sodium  hydroxide  in  the  displacing 
fluid. 

The  adsorption  phenomenon,  both  experimentally  and 
theoret i ca 1 ly . 
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4.  EXPERIMENTAL  EQUIPMENT,  MATERIALS  AND  PROCEDURE 


4.1  Experimental  Equipment 

The  cone  flooding  apparatus  is  shown  in  Figure  2.  The 
displacement  pump  was  a  constant  injection  rate  Ruska  pump 
with  mercury  as  the  displacing  fluid.  All  lines  connecting 
the  pumps  with  the  stainless  steel  cylinders  were  3.175x10~3 
metre  (1/8  inch)  stainless  steel  tubing.  All  lines  upstream 
of  the  cylinders  were  6.35x10-3  metre  (1/4  inch)  stainless 
steel  tubing. 

The  cylinders  had  capacities  of  approximately  3.6x10-3 
cubic  metres  (3.6  litres)  and  were  fitted  with  pressure 
rupture  discs  at  the  upstream  end.  The  system  was  designed 
for  a  maximum  pressure  of  13.1  MPa  (1900  psi ) . 

A  mercury  manometer  and  a  Heise  pressure  gauge  were 
used  to  measure  the  pressure  drop  across  the  core  pack. 
Pressure  transducers  connected  to  a  chart  recorder  were  also 
used  to  provide  a  pressure  history  for  each  run. 


4.2  Porous  Media 

4.2.1  Sand 

Unconsolidated  sandpacks  were  used  throughout  this 
study.  The  sandpacks  were  an  Ottawa  silica  sand  with  a 
reported  grain  size  range  of  80-120  mesh  (Fisher  Scientific 
S -  1 5 1 ) .  The  sand  was  used  as  received  from  the  supplier. 
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SCHEMATIC  OF  CORE  FLOODING  APPARATUS 
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4.2.2  Coreholders 

The  sandpacKs  were  housed  in  stainless  steel 
coreholders  of  5.04  cm  nominal  diameter.  The  lengths  were 
all  in  the  range  of  one  metre. 

4.2.3  Packing  Procedure 

A  wet  packing  procedure  was  used  in  an  attempt  to 
produce  similar  core  properties  from  run  to  run.  This 
procedure  involved  the  following  steps: 

1)  The  steel  coreholder  was  placed  in  a  vertical 
position  with  the  downstream  endcap  in  place. 

2)  Approximately  seven  centimetres  of  distilled 
water  were  placed  in  the  coreholder.  Sand  and 
water  were  alternately  added  in  an  attempt  to 
insure  the  sand  always  fell  through  the  same 
head  of  water. 

3)  During  the  filling  process  the  coreholder  was 
vibrated  by  an  air  vibrator. 

4)  Addition  continued  until  the  coreholder  was 
full.  The  system  was  then  vibrated  overnight. 

5)  Finally  the  upstream  endcap  was  screwed  in 
place . 

The  core  was  then  ready  for  use. 
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4.3  Fluids 

The  following  fluids  were  used  throughout  the 
exper iments : 

4.3.1  Brine 

An  artificial  brine  was  prepared  using  88.84  Kg/m3  of 
reagent  grade  sodium  chloride  dissolved  in  distilled  water. 
This  was  the  same  brine  as  used  by  Gard i ner ( 1 977 ) , 

Sco 1 1 ( 1965)  and  Scot t ( 1971 ) .  This  brine  was  used  as  both 
connate  water  and  as  the  displacing  fluid  for  the  brine 
f  1  ood . 

4.3.2  Caustic  Solution 

In  light  of  previous  work  done  in  the  area  of  optimum 
caustic  concentration,  0.1  %  by  weight  sodium  hydroxide  in 
brine  was  used  for  the  initial  runs.  This  concentration  was 
judged  to  be  optimal  by  Flock,  et  a  1 . ( 1977)  and 
Farouq  Ali,  et  a  1 . ( 1979) .  For  Runs  8,10,11  and  12  the 
caustic  concentration  was  reduced  to  0.01  %  by  weight  sodium 
hydroxide  in  brine,  and  for  Run  9  was  increased  to  2.0  %  by 
weight.  These  runs  were  performed  in  order  to  invesigate  the 
sensitivity  of  recovery  to  caustic  concentration.  The 
caustic  solutions  were  made  up  from  reagent  grade  sodium 
hydroxide  and  the  brine  solution  described  above. 

4.3.3  Oil 

The  oil  used  in  this  study  was  from  the  Wainwright  B 
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Pool.  Excess  water  was  removed  from  the  received  sample  by 
the  use  of  a  rotary  dryer . 


4.4  Fluid  Properties 

4.4.1  Surface  and  Interfacial  Tensions 

The  interfacial  tensions  between  the  crude  oil  and  the 
various  concentrations  of  sodium  hydroxide  in  brine  were 
obtained  using  the  spinning  drop  apparatus . The  densities 
required  for  the  calculation  of  the  interfacial  tensions 
were  measured  with  a  density  meter. 

The  surface  tensions  were  obtained  with  a  du  Nouy 
tensiometer.  These  measurements  were  performed  at  23° C. 

4.4.2  Viscosities 

The  viscosities  of  the  various  fluids  were  obtained  by 
the  use  of  Canon-Fenske  viscometers.  The  determinations  were 
made  at  23°C. 

Appendix  A  details  the  properties  of  the  various 
fluids. 


4.5  Sandpack  Property  Determination 

Prior  to  saturation  of  the  sandpacks  with  oil,  various 
physical  properties  were  measured.  These  measurements  were 
performed  by  routine  methods  and  consequently  are  only 
briefly  described  below. 
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4.5.1  Poros i ty 

The  porosity  of  the  pack  was  determined  via  a  material 
balance  approach  using  the  displacement  of  distilled  water 
by  brine. 

Once  the  pack  was  fully  saturated  with  distilled  water, 
the  water  was  displaced  horizontally,  at  a  constant  rate,  by 
brine,  and  the  chloride  content  of  the  effluent  monitored  by 
refractive  index.  Once  the  chloride  content  of  the  effluent 
reached  that  of  the  displacing  fluid,  the  displacement  was 
stopped.  A  plot,  such  as  Figure  3,  was  constructed  and  the 
area  above  the  curve,  obtained  with  a  planimeter,  yielded 
the  displaced  volume.  The  pore  volume  was  then  evaluated  by 
subtraction  of  the  endcap  volumes. 

4.5.2  Permeabi 1 i ty 

The  absolute  permeability  was  obtained  by  flowing 
distilled  water  through  the  horizontal  pack  at  various  rates 
and  utilizing  Darcy's  Linear  Flow  Equation. 

The  physical  properties  of  the  cores  are  provided  in 
Appendix  B. 

4.5.3  Fluid  Saturations 

Once  the  sandpack  had  been  saturated  with  the  brine 
solution  it  was  suspended  vertically.  Crude  oil  was  then 
injected  down  through  the  pack  and  the  effluent  collected. 
When  the  effluent  contained  less  than  1%  water  the  pack  was 
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assumed  to  be  at  the  irreducible  water  saturation.  The  fluid 
saturations  were  then  calculated  via  material  balance. 


4.6  Displacement  Procedure 

The  displacement  tests  were  run  at  flow  rates  from 
6.94x10-9  m3/sec  to  3.33x10-6  m3/sec  (2.5  to  1200  cc/hour ) , 
which  correspond  to  apparent  linear  rates  of  between  0.003 
m/day  and  14.4  m/day  (0.01  to  47.6  feet/day) .  The  produced 
effluent  was  collected  in  50  ml  centrifuge  tubes  by  an 
indexing  sample  collector.  Where  required,  the  water  and  oil 
phases  were  separated  by  the  addition  of  toluene  and 
subsequent  centrifuging.  The  pH  of  the  aqueous  phase  was 
determined  by  a  digital  pH  meter.  Table  1  summarizes  the 
initial  conditions  of  the  cores  prior  to  the  commencement  of 
the  runs.  Table  2  provides  a  summary  of  the  conditions  under 
which  the  runs  were  performed. 

In  an  attempt  to  reduce  the  effect  of  gravity  in  the 
very  slow  floods,  Runs  4,6  and  11  were  performed  vertically. 
The  criterion  that  only  rates  less  than  0.3  m/day  need  be 
done  vertically  was  arbitrary. 

A  brine  flood  (Run  1)  was  performed  as  a  standard  for 
the  sodium  hydroxide  floods.  This  run  was  done  at  an 
apparent  linear  velocity  of  0.48  m/day  and  the  complete 
results  are  presented  in  Appendix  C. 

Run  3  involved  the  use  of  a  "soak  period"  of  100  hours. 
This  was  done  by  first  injecting  0.25  pore  volumes  of  0.1  % 
by  weight  sodium  hydroxide  in  brine.  The  core  was  then 
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Table  1 


Initial  Conditions  of  Cores 


Ini t i al  Ini t i al 


Run  # 

Core  # 

Di spl aci ng 
Fluid 

Water  in 
(ml  ) 

P  1  ace 

(%pv) 

Oi  1  in 
(ml  ) 

P  1  ace 

( %pv) 

1 

3 

Br  i  ne 

63.2 

8.3 

698.5 

91.7 

2 

3 

Caustic 

64.7 

8.5 

697.0 

91  .5 

3 

2 

Caus  t i c 

108.1 

15.4 

593 . 6 

84.6 

a 

1 

Caustic 

115.0 

14.6 

675.0 

85.4 

5 

4 

Caustic 

94.0 

13.0 

630.0 

87.0 

6 

3 

Caustic 

117.9 

14.5 

693.0 

85.5 

7 

2 

Caustic 

76.9 

11.1 

618.6 

88.9 

8 

4 

Caustic 

72.9 

10.1 

649.7 

89.9 

9 

3 

Caust i c 

94.5 

11.7 

715.4 

88.3 

10 

1 

Caust i c 

81  .3 

11.3 

637.  1 

88.7 

1 1 

4 

Caustic 

74.  1 

10.5 

632.4 

89.5 

12 

2 

Caust ic 

91  .7 

12.4 

650.5 

87.6 
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Table  2 


Run  Summary 


Volumetric  Apparent  Injected  Caustic 

Flow  Rate  Linear  Rate  Concentration 


Run  0 

( cc/hour ) 

(m3/sec ) 

(m/day ) 

( feet/day ) 

(Wt  % 

1 

40.0 

1  .  1  x  1  0  - 7 

0.48 

1.58 

0.00 

2 

40.0 

1  .  1  x  1  0  - 7 

0.48 

1  .58 

0.  10 

*  3 

40.0 

1  .  1  x  1  0  - 7 

0.48 

1.58 

0.  10 

0  4 

2.5 

6 . 9x 1 0 " 9 

0.03 

0.  10 

0.  10 

5 

1200.0 

3.3x10"® 

14.49 

47.55 

0.10 

0  6 

12.5 

3.5x10‘8 

0.  15 

0.49 

0.  10 

7 

250.0 

6 . 9x 1 0 " 7 

3.01 

9.87 

0.10 

8 

250.0 

6 . 9x  1  0 " 7 

3.02 

9.91 

0.01 

9 

250.0 

6 . 9x 1 0" 7 

3.01 

9.87 

2.00 

10 

40.0 

1 . IxlO-7 

0.48 

1.58 

0.01 

0  1  1 

12.5 

3.5x10-® 

0.15 

0.49 

0.01 

12 

1200.0 

3.3x10-® 

14.49 

47.55 

0.01 

+  13 

1200.0 

3.3x10-® 

14.49 

47.55 

0.10 

*  This  run  involved  a  soak  period. 

0  These  floods  were  performed  vertically  upwards. 
+  This  was  the  caustic  displacement  of  brine. 
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sealed  and  allowed  to  soak  for  100  hours  without  rotation. 
After  this  period  the  displacement  was  restarted  and 
continued  until  2  pore  volumes  had  been  injected. 

Run  13  was  designed  to  study  the  adsorption  phenomenon, 
and  involved  the  displacement  of  brine  by  0 . 1  %  by  weight 
sodium  hydroxide  in  brine.  The  core  was  prepared  in  the 
normal  way  and  the  displacement  of  brine  was  carried  out  at 
3.3x10*6  m3/sec  (1200  cc/hour). 


5.  DISCUSSION  OF  RESULTS 

The  experimental  part  of  this  study  involved  thirteen  runs. 
The  theoretical  work  involved  the  determination  of  the 
adsorption  kinetics  parameters,  and  a  predicted 
concentration  profile.  Run  1  was  a  brine  flood  for  use  as  as 
a  reference  run  for  the  caustic  tests.  Runs  2  to  7  involved 
0.1  %  by  weight  caustic  displacing  the  Wainwright  crude  at 
various  different  linear  rates.  Runs  8,  and  10  to  12  were 
performed  with  0.01  %  by  weight  sodium  hydroxide  in  brine. 
Run  13  was  a  simple  displacement  of  brine  by  0 . 1  %  by  weight 
sodium  hydroxide  in  brine,  with  a  view  to  testing  the 
adsorption  theory. 


5.1  Sandpack  Properties 

The  wet  packing  procedure  employed  during  this  study, 
discussed  in  Section  4.2.3,  was  expected  to  provide 
reproduci bi 1 i ty  in  sandpack  properties  such  as,  porosity  and 
absolute  permeability.  The  average  porosity  obtained  was 
34.6  %  with  a  standard  deviation  of  1.2  %.  The  average 
absolute  permeability  was  19.3  darcys  with  a  standard 
deviation  of  1.7.  A  summary  of  the  physical  properties  of 
the  individual  sandpacks  is  provided  in  Table  B-1. 

The  undesireably  large  variations  observed  in  the 
physical  properties  reflect  the  need  for  standardizing  the 
packing  procedure.  More  important,  one  should  insure  that 
all  packs  are  made,  not  only  with  the  same  grade  of  sand 
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but,  from  the  same  manufacturers'  lot  number.  This  point  is 
emphasized  by  the  fact  that  the  sandpacKs  for  Runs  4  and  6, 
which  yielded  high  values  of  porosity  and  permeability, 
employed  sand  from  a  different  lot  than  the  one  used  for 
most  of  the  other  runs.  The  sandpacK  for  Run  13,  which 
yielded  a  below  average  permeability  was  from  a  third  lot. 


5.2  Initial  Fluid  Saturations 

The  initial  saturations  of  the  sandpacKs  are  detailed 
in  Table  1  in  the  previous  chapter.  The  average  initial 
water  saturation  obtained  in  the  sandpacKs  for  Runs  1  to  12 
was  11.8  %  with  a  standard  deviation  of  2.2  %. 

The  variation  in  the  initial  water  saturation  was 
greater  than  expected,  and  is  again  partly  attributed  to  the 
variations  in  the  sand  in  the  different  lots.  The  sandpacKs 
for  Runs  4  and  6,  for  example,  had  initial  water  saturations 
of  14.6  %  and  14.5  %,  respectively,  2.8  %  and  2.7  %  greater 
than  the  average.  The  sandpacK  for  Run  3,  however,  exhibited 
another  high  value  despite  the  fact  that  the  porosity  and 
permeability  were  close  to  the  average  values. 

In  all,  however,  the  fluid  saturations  are  acceptable 
for  quantitative  comparison  of  the  runs.  The  normalizing  of 
the  water  saturation  in  the  relative  permeability  ratio 
calculations  is  an  attempt  to  further  reduce  the  effect  of 


these  variations. 
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5.3  Interfacial  Tension 

The  interfacial  tension  data  between  various  sodium 
hydroxide  solutions  and  the  Wainwright  crude  are  shown  in 
Figure  4.  It  is  evident  from  Figure  4  that  there  is  no  well 
defined  minimum  in  interfacial  tension  in  the  range  of 
caustic  concentrations  studied. 

The  interfacial  tension  is,  however,  reduced  by  the 
addition  of  sodium  hydroxide.  There  is  a  low  value  of 
interfacial  tension  at  0.1  %  by  weight  and  this  was  taken  as 
sufficient  evidence  to  use  this  concentration  for  the 
initial  tests . 

This  finding  is  in  conflict  with  the  work  of 
Chan(1979),  who  found  a  distinct  minimum  in  the  interfacial 
tension  at  0.1  %  by  weight.  Although  Chan(1979)  used  the 
Wainwright  crude,  it  should  be  noted  that  the  viscosity  was 
markedly  different,  800  mPa.s  versus  408  mPa.s.  He  also 
employed  formation  water,  while  this  study  employed  an 
artificial  brine. 


5.4  Displacement  Tests 

A  summary  of  the  displacement  tests  has  been  given  as 
Table  2  in  Chapter  4.  For  the  purpose  of  initial  discussion 
the  displacement  tests  have  been  divided  into  five  sections: 

1)  Run  1 :  A  brine  flood  to  be  used  as  a  base  run 
for  comparing  the  caustic  floods. 

2)  Runs  2-7:  These  runs  were  performed  using 


0.1  %  by  weight  sodium  hydroxide  in  brine,  in 
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order  to  investigate  the  effect  of  rate  on  oil 
recovery.  Run  3  also  studied  the  effect  of  a 
soak  period. 

3)  Runs  8,  10  -  12:  These  runs  investigated  the 
rate  effect  for  a  different  caustic 
concentration,  0.01  %  by  weight  sodium  hydroxide 
in  brine. 

4)  Run  9 :  This  run  was  for  a  sodium  hydroxide 
concentration  of  2.0  %  by  weight.  This  allowed 
comparison  of  the  oil  recovery  at  a  particular 
rate  for  3  different  caustic  concentrations. 

5)  Run  1 3 :  This  run  involved  the  displacement  of 
brine  by  0 . 1  %  by  weight  sodium  hydroxide  in 
brine,  to  obtain  an  experimental  caustic 
concentration  profile  in  the  effluent  for 
comparison  with  that  predicted  by  the  theory. 

5.4.1  Brine  F lood 

Run  1  involved  the  flooding  of  a  core  with  an 
artificial  brine  containing  88.84  kg/m3  sodium  chloride.  The 
production  history  for  this  run  is  provided  as  Table  C-1  and 
Figure  C-1  in  Appendix  C.  This  test  was  performed  at  an 
apparent  linear  rate  of  0.48  m/day.  The  recovery  at  2  pore 
volumes  injected  was  53.9  %  of  the  initial  oil  in  place. 

The  behaviour  of  the  flood  was  as  classically  expected, 
and  the  plot  of  relative  permeability  ratio  versus  a 
normalized  water  saturation  produced  the  expected  straight 
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line  relationship.  This  graph  is  presented  as  Figure  03  in 
Appendix  C. 

The  calculation  of  the  relative  permeability  ratio  was 
performed  using  the  Welge  integration  of  the 
Buckley-Leveret t  equation  from  displacement  data,  as  given 
by  Co  1 1 i ns ( 1 96 1 ) .  The  decision  to  use  a  normalized  water 
saturation  rather  than  the  actual  saturation,  as  suggested 
by  Col  1 i ns ( 1 96 1 ) ,  was  based  on  the  desire  to  quantitatively 
compare  the  relative  permeability  ratio,  from  the  various 
tests,  on  a  consistent  basis.  Due  to  the  lack  of  knowledge 
of  the  residual  oil  saturation,  the  following  definition  of 
normalized  saturation  was  used: 


wi 


5.4.2  Tests  with  0.1  %  by  Weight  Sodium  Hydroxide 

The  initial  caustic  floods  were  performed  using  0.1  % 
by  weight  sodium  hydroxide  in  brine  as  the  displacing  fluid. 
Runs  2  and  3  were  performed  at  the  same  linear  rate 
(0.48  m/day)  and  differed  only  in  that  Run  3  utilized  a 
" soak  per i od" . 

Run  3  was  performed  to  examine  whether  the  "real  time" 
concept  proposed  by  Chan ( 1979)  plays  a  role  in  the  system 
under  study.  It  is  apparent  from  Figure  5  that  the  recovery 
histories  for  the  two  runs  were  essentially  the  same. 

Chan ( 1979)  found  that  the  recovery  improved  dramatically, 
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once  the  contact  time  exceeded  some  minimum.  As  this  did  not 
appear  to  be  the  case  for  the  system  under  study  the  concept 
of  a  time  dependence  was  not  studied  further. 

The  remaining  0.1  %  by  weight  floods,  Runs  4  to  7, 
studied  the  effect  of  rate  on  oil  recovery  at  both 
breakthrough,  and  two  pore  volumes  injected.  The  complete 
series  of  production  histories  is  shown  on  Figure  6.  The 
individual  run  histories  are  provided  in  both  tabular  and 
graphical  form  in  Appendix  C.  Figure  6  demonstrates  the 
apparent  insensitivity  of  the  recovery  process  to  rate.  All 
the  tests,  with  the  exception  of  Run  5,  follow  an  almost 
identical  path.  The  breakthrough  in  Run  5,  performed  at 
14.49  m/day  (47.55  ft/day) ,  was  earlier  than  in  the  other 
runs,  but  the  recovery  at  2  pore  volumes  injected  is 
difficult  to  distinguish  from  the  other  runs. 

This  data  is  more  elegantly  displayed  in  Figure  7  from 
which  it  may  be  clearly  seen  that  there  is  little 
significant  difference  in  oil  recovery  except  at  the 
extremes  of  linear  velocity. 

This  would  suggest  that  the  dominant  recovery  mechanism 
is  largely  insensitive  to  rate.  This  will  be  discussed 
further  in  Section  5.5  . 

5.4.3  Tests  with  0.01  %  bv  Weight  Sodium  Hydroxide 

These  floods  were  undertaken  to  invesigate  the 
sensitivity  of  the  dominant  displacement  mechanism  to 
caustic  concentration.  Also,  by  performing  a  series  of  tests 
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at  different  apparent  linear  rates,  the  rate  sensitivity 
could  be  further  investigated. 

It  is  immediately  clear  from  Figure  8  that  the  dominant 
recovery  mechanism  is  significantly  affected  by  the  caustic 
concentration.  This  figure  shows  that,  at  an  apparent  linear 
velocity  of  3.01  m/day  (9.87  ft/day) ,  the  0.01  %  by  weight 
solution  is  superior  to  the  2.0  or  the  0.1  %  by  weight 
solutions.  The  improvement  over  the  0.1  %  case  is  only 
evident  after  the  injection  of  approximately  1.05  pore 
volumes.  As  can  be  seen  from  the  oil  cut  histories,  provided 
in  Appendix  C,  the  breakthrough  in  Run  8,  using  0.01  % 
caustic,  was  earlier  but,  the  oil  cut  was  higher  for  the 
remainder  of  the  displacement.  This  trend  persisted  for  all 
the  rates  studied  but  is  most  apparent  in  this  particular 
case . 

Examination  of  Figure  7  shows  that  in  all  cases  the 
breakthrough  recovery  for  the  0.01  %  by  weight  solution  was 
lower  than  for  the  correspondi ng  0.1  %  case.  The  recovery  at 
2  pore  volumes  injected,  on  the  other  hand,  was  higher  in 
all  cases,  except  Run  11  which  had  essentially  the  same 
value  as  its  counterpart,  Run  6. 

A  further  interesting  point  from  Figure  7  is  that  the 
recovery  appears  to  be  more  rate  sensitive  in  the  0.01  %  by 
weight  system,  than  in  the  0.1  %  by  weight  system.  For  the 
former ,  the  recovery  at  2  pore  volumes  injected  goes  through 
a  maximum  at  approximately  3  m/day.  It  appears  then,  that 
the  floods  performed  with  the  0.01  %  by  weight  caustic 
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involved  a  different  recovery  mechanism  from  that  operating 
in  the  0.1  %  case.  Possible  explanations  for  this  behaviour 
will  be  presented  in  Section  5.5  . 

5.4.4  Tests  with  2.00  %  by  Weight  Sodium  Hydroxide 

This  run  was  designed  to  further  examine  the  effect  of 
caustic  concentration  on  the  recovery  from  the 
sand/oi 1 /water  system  under  study.  It  is  apparent  from  both 
Figure  7  and  Figure  8  that  there  was  little  significant 
improvement  in  oil  recovery,  over  the  brine  flood,  when 
using  the  2.0  %  by  weight  solution. 

An  examination  of  the  respective  relative  permeability 
ratio  plots  for  the  brine  and  2.0  %  caustic  floods, 

Figures  03  and  027,  shows  that  there  was  no  major  shift  in 
the  relative  permeability  ratio,  as  was  experienced  with  the 
other  two  caustic  concentrations.  This  information  is 
presented  in  Figure  9.  It  is  also  apparent  from  this  figure 
that  despite  the  fact  that  breakthrough  occurred  later  in 
the  2.0  %  by  weight  case,  there  is  little  or  no  improvement 
in  relative  permeability  ratio  beyond  a  normalized 
saturation  of  approximately  0.3  .  This  is  also  evident  from 
the  recovery  curves  in  Figure  8  where  Runs  1  and  9  provide 

much  the  same  production  history. 

It  appears  then,  from  Run  9,  that  an  excessively  large 
caustic  concentration  is  in  fact  detrimental  to  recovery 
enhancement.  Subkow(1942)  stated  that  high  electolytic 
concentrations  could  prevent  the  formation  of  oil-in-water 


%:  ,  *  + .  3  mwjm  X  ?■ 


RELATIVE  PERMEABILITY  RATIO  (Krw/Kro) 


49 


FIG  9  :  RELATIVE 

o  permeability  RATIO  CURVES 


50 

emulsions,  which  is  an  essential  step  in  both  the 
emulsification  and  entrapment,  and  emulsification  and 
entrainment  recovery  mechanisms. 

Emulsion  formation,  reduction  of  the  interfacial 
tension  at  the  oil-water  interface  and  wettability 
alterations  will  all  cause  a  shift  in  the  relative 
permeability  ratio  towards  a  less  favourable  water  mobility. 
The  fact  that,  for  the  2.0  %  case,  no  shift  occurred 
suggests  that  none  of  these  mechanisms  were  present  in  this 
f 1 ood . 


5.5  Comparison  of  the  Effect  of  Rate  on  the  0.1  and  0.01  % 

by  Weight  Tests 

The  effect  of  rate  on  oil  recovery  may  be  due  to  any, 
or  all,  of  the  following  factors: 

1 )  Emu  1 s i ons 

2)  Viscous  forces 

3 )  Capi 1 1  ary  forces 

4)  Adsorption 

5)  Gravity 

The  displacements  at  low  flow  rates  were  performed 
vertically  in  an  attempt  to  negate  the  influence  of  gravity. 
There  is,  however,  little  that  can  be  done  to  control  any  of 
the  other  factors. 

A  study  of  total  caustic  consumption  in  the  0.1  %  by 
weight  floods  shows  that  the  consumption,  at  this 
concentration,  is  insensitive  to  rate.  This  total  caustic 
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consumption  incorporates  both  adsorption  on  the  sand,  and 
reactions  with  the  oil  and  sand.  There  is  some  scatter  in 
the  data,  which  can  be  seen  on  Figure  10,  but  the  general 
trend  is  obvious.  It  thus  appears  that  the  system  has 
reached  an  equilibrium  with  the  flowing  caustic  stream.  The 
consumption  of  caustic  by  reaction  with  the  oil  appears  to 
be  the  major  element  in  the  total  consumption.  The 
consumption  in  Run  13,  shown  as  *  on  Figure  10,  when  no  oil 
was  present,  is  a  full  order  of  magnitude  less  than  the 
correspondi ng  run  with  oil  present. 

It  was  noted  that  the  apparent  linear  rate  had  markedly 
different  effects  in  the  0.01  %  by  weight  and  0.1  %  by 
weight  cases.  The  differing  effects  of  rate  may  be  seen  in 
Figure  11,  where  the  character  of  the  two  curves  is  notably 
different,  with  the  0.1  %  by  weight  flood  breaking  through 
later  and  the  percentage  of  oil  in  the  effluent  dropping  off 
quicker  than  in  the  0.01  %  by  weight  case.  The  recoveries  at 
two  pore  volumes  injected  are  similiar,  61.4  %  for  the  0.1 
and  61.2  %  for  the  0.01  %  by  weight  case.  This,  despite  the 
markedly  different  oil  cut  histories.  This  fact  suggests 
different  recovery  mechanisms  are  in  play  for  the  two  cases. 

An  examination  of  the  relative  permeability  ratio 
curves  on  Figure  9  shows  that  the  two  cases  produce  curves 
of  different  character.  The  0.01  %  by  weight  displacement 
produces  an  "S"  shaped  curve,  while  the  0.1  %  by  weight 
displacement  yields  a  slighty  concave  curve.  In  both  cases, 
when  compared  with  the  brine  curve,  the  curves  are  shifted 
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towards  a  less  favourable  water  mobility.  The  crossover  of 
the  two  caustic  curves  shows  that  early  in  the  life  of  the 
displacement  the  0.1  %  case  is  superior;  but  in  the  latter 
stages,  the  0.01  %  case  is  superior. 

A  possible  explanation  for  this  behaviour  may  be  found 
in  the  mechanisms  in  play.  It  appears  that  in  the  0.1  % 
system,  the  improved  oil  recovery  is  due  to  emulsification 
and  entrainment.  The  interfacial  tension  is  low  in  this 
system,  and  the  oil  will  be  emulsified  in  situ  and  carried 
with  the  flowing  alkali  and  then  produced.  The  fact  that  no 
emulsions  were  observed  in  the  effluent  does  not  rule  out 
this  mechanism.  The  emulsions  may  have  been  so  unstable  that 
they  broke  down  in  the  collection  tubes.  It  was  noted  during 
the  displacements  that  the  effluent  took  a  finite  time  to 
separate  into  distinct  oil  and  water  phases. 

In  the  0.01  %  case,  on  the  other  hand,  the  dominant 
mechanism  appears  to  be  emulsification  and  entrapment.  The 
interfacial  tension  in  this  system  was  higher  than  that  in 
the  0.1  %  case  and,  as  a  result,  the  emulsified  oil  could 
not  be  pushed  through  the  pore  throats.  It  appears  from  the 
relative  permeability  ratio  curve  on  Figure  9  that,  the  oil 
droplets  block  off  the  swept  region,  causing  the  decreased 
water  mobility  demonstrated  by  the  flattened  section  of  the 
curve  at  intermediate  values  of  normalized  saturation.  This 
results  in  an  increased  displacement  efficiency.  The 
emulsified  oil  in  the  entrainment  mechanism  does  not  block 
off  the  pore  throats  and,  as  a  result,  the  oil  cut  drops  off 
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rapidly  after  breakthrough . 

The  differing  effects  of  linear  rate  on  recovery  are 
due  to  these  different  mechanisms.  The  entrainment  mechanism 
is  essentially  independent  of  rate,  as  the  entrainment 
process  itself  is  rate  independent.  The  differences  in 
recovery,  at  the  extreme  rates,  are  simply  due  to  the 
differences  in  displacement  efficiency,  which  are  caused  by 
the  viscous  and  capillary  forces.  In  the  entrapment 
mechanism,  however,  the  linear  rate  effects  the  recovery  in 
two  ways : 

1)  The  increased  rates,  under  increased  gradients, 
will  push  the  oil  droplets  through  the  pore 
throats,  thereby  preventing  the  blocking  of 
these  throats. 

2)  At  reduced  rates,  the  caustic  consumption  will 
result  in  the  interfacial  tension  not  being 
sufficiently  reduced.  This  factor  is  not  crucial 
in  the  entrainment  mechanism  because  of  the 
magnitude  of  the  input  concentration. 

In  the  case  of  the  entrapment  mechanism,  less  sodium 
hydroxide  should  be  produced  as  the  displacing  fluid 
contacts  more  of  the  sandpack.  The  pH  of  the  effluent  was 
monitored  and  it  is  evident  that  the  majority  of  the  caustic 
was  being  consumed.  Examination  of  experimental  results,  in 
Tables  C-1  to  C-12,  shows  a  maximum  value  for  C/CQ  of  0.285 
(Run  10).  This  compares  with  a  minimum  value  of  0.530  in  the 
entrainment  displacements. 
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5.6  Adsorption  Study 

Run  13  was  undertaken  to  quantitatively  study  the 
phenomena  of  caustic  dispersion  and  adsorption  in  the 
sandpack.  The  flood  was  conducted  using  0.1  %  by  weight 
sodium  hydroxide  in  brine  solution  displacing  brine.  Prior 
to  the  caustic  flood,  the  displacement  of  distilled  water  by 
brine,  used  in  the  porosity  determination,  was  used  to 
evaluate  the  dispersion  coefficient,  D,  by  the  method 
proposed  by  Brigham,  et  a  1 . ( 1 96 1 ) . 

As  a  check  on  the  value  of  D,  the  profile  predicted  by 
Equation  1  was  plotted  on  Figure  12  with  the  experimental 
data.  It  is  apparent  from  this  figure  that  the  predicted 
profile  is  a  resonable  fit  up  to  high  values  of  normalized 
concentration.  The  experimental  data  exhibits  a  "tailing 
off"  not  predicted  by  the  dispersion  equation.  This  "tail" 
is  due  to  dead-end  spaces  in  the  sandpack,  from  which  the 
distilled  water  will  diffuse  as  the  brine  bypasses  it. 
Therefore  it  is  not  pushed  ahead  by  the  brine  as  assumed  by 
the  theory.  The  experimental  and  theoretical  results  are 
also  presented  in  tabular  form  in  Appendix  D. 

Following  the  caustic  flood  the  rate  constant  K-j  was 
obtained  by  point  matching  as  described  in  Section  2. 5. 4. 2  . 
This  procedure  to  obtain  K,  has  a  serious  shortcoming.  In 
this  system  the  value  of  C/CQ  at  50  =1.0  is  very  small,  the 
experimental  value  being  less  than  0.02  .  At  these  very  low 
concentrations  the  accuracy  of  the  readings  would  be 
questionable  and,  as  a  result,  the  value  of  K-j  must  also  be 
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in  doubt.  This  problem  is  eased  somewhat  by  the  relationship 
in  Equation  32  which  relates  the  value  of  K-j  to  the  natural 
log  of  the  concentration.  A  functional  relationship  such  as 
this  tends  to  smooth  out  possible  errors. 

The  experimental  caustic  concentration  profile  is  shown 
in  Figure  13.  In  the  case  of  zero  adsorption  the  profile 
would  be  identical  to  that  obtained  in  the  brine 
displacement  of  distilled  water.  Comparison  of  Figures  12 
and  13  shows  this  is  not  the  case.  The  concentration  profile 
has  been  shifted  to  the  right.  This  delay  in  the  caustic 
breakthrough  has  been  caused  by  adsorption  of  the  caustic 
onto  the  rock  surfaces. 

The  next  step  was  to  attempt  to  fit  this  experimental 

profile  to  the  theoretical  models.  An  attempt  was  first  made 

to  use  the  Case  I  solution  of  Banks  and  A 1 i ( 1 964 ) ,  given  as 

Equations  (23)  to  (26).  The  resulting  prediction  is  shown  in 

Figure  13.  It  appears  from  this  figure  that  the  prediction 

is  resonable  for  Co  <  1.0  .  This  finding  is  in  agreement 

with  that  of  Banks  and  A 1 i ( 1 964 ) .  A  possible  explanation  for 

the  leveling  of  the  predicted  profile  is  that,  with  ^  =  0.0 

in  the  rate  equation,  there  is  no  reduction  in  adsorption 

rate  as  C  increases.  The  result  is,  then,  that  C  will 
s  b 

increase  with  time  until  some  equilibrium  is  reached  with 
the  passing  stream.  This  value  is,  however ,  much  higher  than 
that  obtained  in  reality. 

Following  the  failure  of  the  Case  I  model  to  accurately 
predict  the  experimental  results,  the  logical  place  to  turn 
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was  the  general  solution,  Equation  (20).  A  value  for  K2  was 
obtained  by  point  matching,  but  attempts  to  use  the  general 
solution  failed  due  to  mathematical  difficulties.  It  was 
found  that  the  exponential  terms  of  Equations  (20)  and  (21) 
generated  numbers  of  such  extreme  magnitude  that  computation 
became  impossible.  Attempts  were  made  to  rewrite  the 
equations  in  more  manageable  forms,  but  these  met  with 
little  success . 

Due  to  the  obvious  effects  of  dispersion,  the  Case  II 
model  of  Banks  and  A 1 i ( 1 964 )  was  not  used.  This  model 
assumes  that  there  is  no  dispersion.  It  is  proposed  that  the 
di spersion-adsorpt ion  equation  be  solved  numerically,  using 
a  more  appropriate  adsorption  relationship. 


5.7  Experimental  Accuracy 

It  is  important  in  quantitative  research,  such  as  this 
study,  to  fully  appreciate  the  limitations  of  the 
experimental  results.  The  data  collected  in  this  study  may 
be  divided  into  2  sections: 

1)  Volumetric  measurements  of  the  effluent. 

2)  Measurement  of  the  effluent  concentrations. 

The  effluent  from  the  displacement  tests  was  collected 
in  50  ml  centrifuge  tubes.  These  tubes  allow  accuracy  of  no 
better  than  ±  0.5  ml.  In  the  early  stages  of  the 
displacements,  with  large  oil  volumes,  this  represents  a 
possible  error  of  1%.  In  the  later  stages,  however,  when  oil 
volumes  are  reduced  to  as  little  as  2/4  of  the  effluent,  this 
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possible  error  is  much  more  dramatic.  In  the  case  of  the 
production  history,  any  error  is  damped  out  by  the  large 
cumulative  production  but,  in  the  calculation  of  the 
relative  permeability  ratio,  it  is  not.  This  is  due  to  these 
calculations  using  the  slope  of  the  production  history 
between  individual  points.  To  overcome  this,  these 
calculations  were  performed  on  various  average  points, 
obtained  by  fitting  a  straight  line  to  5  data  points.  The 
straight  line  sections  were  overlapped  to  provide 
cont i nui ty . 

The  brine  concentration  in  the  effluent  (from  the 
porosity  determination  displacements)  was  monitored  by 
refractive  index.  The  Spencer  Ref ractometer  used,  allowed 
readings  to  ±  0.0005.  This  represents  a  possible  error  of 
±  3.4%  in  the  resulting  concentrations. 

The  caustic  concentrations  could  not  be  monitored  by 
refractive  index  as  the  difference  between  the  indices  of 
the  caustic  solutions  and  that  of  the  brine  was  too  small. 
These  concentrations  were  obtained  through  pH.  The 
relationship  between  concentration  and  pH,  as  shown  in 
Figure  14,  is  such  that  the  accuracy  would  vary  from  0.5  % 
to  11  %.  At  low  values  of  sodium  hydroxide  concentration  the 
pH  increases  rapidly  for  small  increases  in  concentration, 
while  at  higher  values  the  pH  increases  much  more  slowly. 
This  means  that  the  possible  error  increases  with  caustic 
concentration.  This  may  lead  to  considerable  errors  in  the 
calculation  of  sodium  hydroxide  consumption. 
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The  consumption  of  sodium  hydroxide  in  Run  13  may  be 
due  to  adsorption  on  and/or  reaction  with  the  silica  sand. 
The  solubility  of  silica  in  hot  alkali  is  reported  by  Perry 
and  Chi  1  ton ( 1 973 )  to  be  "very  slight".  In  the  system  under 
study,  where  the  displacements  were  performed  at  23° C,  the 
solubility,  and  as  a  result  the  caustic  consumption  due  to 
the  silica,  should  be  negligible. 


6.  SUMMARY  AND  CONCLUSIONS 

Laboratory  caustic  floods  were  conducted  in 
unconsolidated  sandpacKs  saturated  with  an  artificial  brine 
and  the  Wainwright  crude.  A  displacement  of  brine  by  caustic 
was  also  performed,  to  examine  the  adsorption  behaviour.  The 
following  conclusions  were  drawn  from  the  experimental 
results  obtained: 

1.  The  addition  of  0.1  and  0.01  %  by  weight  sodium 
hydroxide  to  the  displacing  fluid  improves  oil 
recovery . 

2.  The  dominant  recovery  mechanism  appears  to  be 
governed  by  the  sodium  hydroxide  concentration. 
Emulsification  and  entrainment  is  suggested  as  the 
dominant  recovery  mechanism  in  the  0.1  %  by  weight 
displacements,  while  emulsification  and  entrapment 
is  used  to  explain  the  behaviour  of  the  0.01  %  by 
weight  displacements. 

3.  The  dominant  recovery  mechanism  in  the  0.1  %  by 
weight  displacements  seems  to  be  essentially 
independent  of  rate,  within  the  framework  of  this 

study. 

4.  The  dominant  recovery  mechanism  in  the  0.01  %  by 
weight  displacements  appears  to  be  rate  dependent. 

5.  In  the  0.1  %  by  weight  displacements  the  total 
caustic  consumption  appeared  insensitive  to  rate. 

The  oil  present  consumed  a  greater  proportion  of  the 
total  than  the  sandpack. 
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6.  Consumption  of  caustic,  by  adsorption  on  and/or 
reaction  with  the  rock  did  occur  during  the 
displacement  of  brine  by  0 . 1  %  by  weight  caustic. 


. 


7.  RECOMMENDATIONS 


The  following  possible  areas  of  study  are  proposed  to 
extend  this  study: 

1.  The  magnitude  of  sodium  hydroxide  adsorption  should 
be  examined  in  relation  to  both  the  input 
concentration  and  the  linear  displacement  rate. 

2.  An  attempt  should  be  made  to  improve  the 
mathematical  representation  of  the  adsorption 
phenomenon,  with  special  attention  to  the  rate 
equation . 

3.  The  effect  of  rate  on  smaller  sodium  hydroxide 
concentrations  should  be  investigated. 

4.  The  sensitivity  of  both  the  adsorption  phenomenon 
and  the  rate  effect  to  rock  type  should  be  examined. 
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Table  A-  1 


Properties  of  Aqueous  Fluids 


Fluid 

Dens i ty 
at  23 °  C 
Kg/m3 

V i scos i ty 
at  23 °  C 
mPa .  s 

Surface 

Tension 

N/m 

Br  i  ne 

(88.84  Kg/m3  NaCl 
in  distilled  water . ) 

1017.4 

1.149 

6.  12x10-2 

0.01  wt  %  NaOH 
in  brine 

1063.7 

1  .  169 

6.09x10-2 

0.10  wt  %  NaOH 
in  brine 

1063.7 

1  .  160 

6.05x10-2 

2.00  wt  %  NaOH 
in  brine 

1085.9 

1  .  182 

6.04x10-2 

Table 

A-  2 

Properties  of  the 

Dr i ed  Crude  Oi 1 

Water  Content 

0.65  wt  % 

Surface  Tension  =  2 
Viscosity  at  23°C  =  403 
Dens i ty  at  23 °  C  =941 
Acid  Number  =  0 
Asphaltene  content  =  6 


91x10-2  N/m 
3  mPa . S 

3  Kg/m3 

527  mg  KOH/g 

06  % 


r  •'  d 


' 
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APPENDIX  B 


Physical  Properties  of  Core  Packs 


Run 

# 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1  1 

12 

13 


roro-F*— ' ^coj^Mco-ti- ^rococo 
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Core 

# 


Table  B-  1 


Physical  Properties  of  Core  Packs 


Length 

Bulk 

Volume 

(cm) 

(cc) 

113.0 

2254.4 

113.0 

2254.4 

105.1 

2096.8 

105.1 

2088.5 

105.7 

2100.4 

113.0 

2254.4 

105.1 

2096.8 

105.7 

2100.4 

113.0 

2254.4 

105.1 

2096.8 

105.9 

2104.4 

106.4 

2122.7 

105.  1 

2096.8 

Pore 

Volume 

Poros i ty 

( cc ) 

(%) 

761.7 

33.8 

761  .7 

33.8 

701.7 

33.5 

790.0 

37.8 

724.0 

34.5 

810.9 

36.0 

695.5 

33.2 

722.6 

34.4 

809.9 

35.9 

718.4 

34.3 

706.5 

33.6 

742.2 

35.0 

709.2 

33.6 

Absolute 
Permeabi 1 i ty 
(darcys ) 

17.8 

17.8 
18.  1 

22.9 

18.5 

22.4 

19.9 
18.  1 
20.  1 

19.4 
17.7 
19.0 

16.6 
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Summary  of  Individual  Displacement  Tests 
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FIG  C . 30  :  RELATIVE 
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FIG  C.3S  :  RELATIVE 
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Table  D-  1 


Comparison  of  Measured  and  Predicted 

Normalized  Brine  Concentration  Profiles 


Predicted  Profile  from  Equation  1 
Experimental  Data  from  Run  13 


Dispersion  Coefficient:  2.312x10-7  m2/sec 
Time  to  Inject  1  PV  :  2128.6  sec 

Pore  Volume  :  709.2  ml 


Pore  Volumes 
Injected 

0.890 
0.949 
0.969 
0.992 
1.010 
1 .038 
1.054 
1  .072 
1  .093 
1.114 
1  .  141 
1  .  160 
1  .  178 
1  .  196 
1  .239 
1  .250 


Measured 
Concent rat  ion 

0.000 
0.014 
0.077 
0.239 
0.585 
0.796 
0.873 
0.930 
0.937 
0.944 
0.951 
0.958 
0.965 
0.993 
0.993 
1  .000 


Ca leu  1 ated 
Concentration 

0.000 
0.040 
0.  146 
0.394 
0.631 
0.894 
0.961 
0.990 
0.999 
1.000 
1.000 
1.000 
1.000 
1  .000 
1.000 
1  .000 
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Table  D-  2 


Comparison  of  Measured  and  Predicted 
Normalized  Caustic  Concentration  Profiles 


Predicted  Profile  from  Equation  23 
Experimental  Data  from  Run  13 


Dispersion  Coefficient:  2.312x10'7  m2/sec 
Rate  Constant  K-j  :  9.750x10'4  sec-1 


Pore  Volumes 
Injected 

0.849 
0.908 
0.930 
0.953 
0.973 
0.996 
1  .017 
1  .037 
1  .058 
1  .077 
1  .097 
1  .  121 
1.148 
1.169 
1.193 
1.215 
1  .239 
1  .272 
1  .561 
1.778 


Measured 
Concentrat i on 


0.000 

0.001 

0.001 

0.003 

0.005 

0.007 

0.027 

0.106 

0.205 

0.295 

0.379 

0.464 

0.566 

0.663 

0.783 

0.819 

0.843 

0.892 

0.952 

0.976 


Ca 1 cu 1 ated 
Concentrat i on 


0.000 

0.000 

0.000 

0.001 

0.004 

0.009 

0.014 

0.016 

0.018 

0.018 

0.018 

0.018 

0.018 

0.018 

0.018 

0.018 

0.018 

0.018 

0.018 

0.018 


